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ABSTRACT 

We present high resolution (i? = 24, 000) L-band spectra of the young intermediate 
mass star V1331 Cyg obtained with NIRSPEC on the Keck II telescope. The spectra 
show strong, rich emission from water and OH that likely arises from the warm surface 



regio n of the circumstellar disk. We explore the use of the new BT2 ([Barber et al 



20061 ) water line list in fitting th e spectra, and we fin d that it does a much better job 



than the well-known HITRAN (jRothman et al.l Il998l ) water line list in the observed 



wavelength range and for the warm temperatures probed by our data. By comparing 
the observed spectra with synthetic disk emission models, we find that the water and 
OH emission lines have similar widths (FWHM ~ 18kms~i). If the line widths are set 
by disk rotation, the OH and water emission lines probe a similar range of disk radii in 
this source. The water and OH emission are consistent with thermal emission for both 
components at a temperature ~ 1500 K. The column densities of the emitting water 
and OH are large, ~ lO^^ cm"^ and ~ lO^'' cm"^, respectively. Such a high column 
density of water is more than adequate to shield the disk midplane from external UV 
irradiation in the event of complete dust settling out of the disk atmosphere, enabling 
chemical synthesis to continue in the midplane despite a harsh external UV environment. 
The large OH-to-wate r ratio is similar to expectations for UV irradiated disks (e.g.. 



Bethell &: Berginll2009l ). although the large OH column density is less easily accounted 



for. 
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1. Introduction 

Near-infrared ro-vibrational emission lines of water and OH are probes of the physical and 



Carr et al. 


2004: 


Naiita et al. 


2007) 



from high accreti on rate systems 



Carr et al 



> 10 ^Mf7) yr ranging from low mass objects such as SVS 13 



2OO4I1 and DG Tau dNaiita et all l200ojL_ ;o higher ni ass objects such as VI 331 Cyg 



(INaiita et all 120091 ') and [BKH2005] 08576nr292 (jThi fc BikI 120051 ') . The water emission in these 



systems is accompanied by 2.3 /im CO overtone emission that also arises from the inner disk. In 
the lower mass objects in which water emission is detected, the emission arises from within ~ 0.5 AU. 
In all systems the emission column densities are large (^10^*^-10^^ cm ) and the temperature of 
the emitting gas is high (~1500K). 

These observations complement studies of water emis sion made at inid-infrared wavelengths 



with the Spitzer Space T e lescope and ground-based facilities (ICarr &: Naiita 



2010a 



Salvk et al. 



2011 



2008 



Pontoppidan et al.ll2010bl: 



Sturm et al. 



Knez et al. 



2011 



2008 



Pontoppidan et al 



20071 ) and longer wavelength 



which measure cooler water, 



studies with the Herschel Space Telescope (e.g., 
probably arising at larger disk radii on average. The mid-infrared water emission arises from sys- 
tems spanning a wid e range of accretion rates, including typical T Tauri stars such as AA Tau 
( Carr &: Naiita 2008 , < 10~ ^ Mq yr~^) and high accretion rate T Tauri stars such as DR Tau and 



AS 205A dSalvk et al.ll2008l . > lO"'^ Mq yr'^). Perhaps because the mid-infrared water emission 
can arise under a wide range of conditions, it appears comrnonly i n Spitzer spectra of T Tauri stars 
(|Carr fc Naiitalboill : Isalvk et al.l[201ll : IPontoppidan et al.lboiOal ). 



Ro-vibrational OH emission is also believed to arise from the inner regions of gaseo us disks. OH 



emission in the L-band has been previou sly reported froin V1331 Cyg and SVS 13 (jNaiita et al 



20071 ). and from DR Tau and AS 205 A (jSalvk et al.ll2008l l. All of these sources also show H2O 
emission in their L-band spectra. OH emission has also been de tected from Herbig stars; the L- 



band spectra of these sources is notably lacking in water emission (jMandell et al.ll2008l : iFedele et al 



20 111 ). The L-band OH emission in the Mandell sources has temperatures of 650-1000 K, and 
column densities of 10^^-10^'' cm~'^. 



^Data presented herein were obtained at the W.M. Keck Observatory from telescope time allocated to the National 
Aeronautics and Space Administration through the agency's scientific partnership with the California Institute of 
Technology and the University of California. The Observatory was made possible by the generous financial support 
of the W.M. Keck Foundation. 
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Here we present high-resolution L-band spectra of V1331 Cyg (LkHa 120), a young star in the 
L988 dark cloud complex. Previous studies of this so urce have suggested a spectral type (A8~G5, 
Kuhilll964l : I Chavarrial 1 1 98 ll : iHamann &: Perssonlll992l ) that is intermediate between those of Herbig 
Ae stars and t ypical T Tauri stars. T here are various estimates of the distance to the source (for 



Herbig &: Dahmll2006l ). with extinction studies of the region toward the L988 dark 



a summary see 

clouds suggestin g a distance of ~ 55 pc. The accretion rate of V1331 Cyg is high. The Br7 line 



flux measured by Eisner et al.l (|2007l ) corresponds to an accretion luminosity of 8.5Lq at a distance 
of 550 pc, putting the accretion rate at > 3 x 10~^-/VfQyr~^. 

As noted above, OH emissio n from V1331 Cyg has been previously reported, based on a small 
portion of the L-band spectrum (jNaiita et al.ll2007l ). The data presented here cover a much larger 
wavelength range, allowing us to constrain the properties of the molecular emission in this spectral 
region. This study complements our earlier study of CO overton e and water emission in the /C-band 
spectrum of V1331 Cyg (see also lNaiita et ahlbood : ICarilll989l l. 



In o ur previous studies of water emission, we used a water line list that was constructed 
from the iPartridge &: Schwenkd (|l997l ) and HITRAN lin e lists and refined empirically in the small 
portion of the i^-band that we stu died (jCarr et al.ll2004l ). Here we explore the use of the new BT2 
water line list (jBarber et al.ll2006l ) in modeling a much larger wavelength region. 



2. Observations 
2.1. Spectroscopy 

High resolution L-band spectra of V1331 Cyg were taken on 11 July 2001 with NIRSPEC 



(jMcLean et al.lll998l ). the multi-order cryogenic echelle facility spectrograph on Keck II. The data 
were obtained by one of us (JRG) through time granted to our NASA/Keck proposal. The spectra 
cover most of the L-band at a resolution of R = X/SX = 24,000 (12.5kms~^) using the 0.432" 
(3-pixel wide) slit. The echelle orders overfill the detector in the dispersion direction, and therefore 
wavelength coverage of most of the L-band was achieved in two grating setups. In the first setup 
(hereafter Ll), the echelle and cross disperser angles were oriented at 62.2° and 33.15°, respectively. 
In the second setup (hereafter L2), the echelle and cross disperser angles were oriented at 64.3° and 
33.29°, respectively. Both grating setups imaged separate portions of orders 20-25 onto the 1024 
X 1024 InSb detector array through the KL filter. Poor telluric transmission in orders 24 and 25 
(at 3.2 /xm and 3.1 iim) prevented our using any data at these grating settings, thus our analysis 
used the spectra we obtained in orders 23-20 (3.27-3.88 /im, Fig. 1). 

To obtain both object and sky spectra, we nodded the telescope along the 24" long slit in 
three positions in an ABC nod pattern, with one position centered on the slit and the two others 
offset by —8" and +8" relative to the slit center. Single exposures of V1331 Cyg were 30 seconds 
X 4 co-adds at each nod position, resulting in 12 minutes of total integration time in ~ 0'.'5 seeing 
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(FWHM L-band) within an airmass range of 1.35-1.45. During the observations, the instrument 
image rotator was kept in stationary mode, keeping the sht physically stationary and thus allowing 
it to slowly rotate on the sky as the alt-az mounted Keck II telescope tracked. 

To remove telluric absorption features, spectra of an early type star (HR 7610, Al IV) were 
taken immediately following the V1331 Cyg observations using the same two grating configurations 
successively over a similar airmass range (1.3-1.4). Spectra of the NIRSPEC internal continuum 
lamp were used for flat-fielding. We used telluric absorption lines identified in the spectra of our 
target and calibration stars for wavelength calibration. 



2.2. Data Reduction 



Initially, bad pixels present in the raw object, telluric, and calibration frames were identified 
and removed by interpolation using the fixpix algorithm within the REDSPEc|^ reduction package 
(a custom echelle reduction package developed by Prato, Kim, and McLean). 



Standard IRAF packages (jMassev et al.l 11992 . 119971 ) were used to reduce the data. Pairs of 
exposures taken at different nod positions closest in time were differenced then divided by the 
normalized internal continuum lamp frame for adequate sky subtraction and flatfielding. Each 
echelle order on the array was then rectified using selected bright sky emission lines (typically 9-12 
lines) present in the object frame, summed over all nod positions. Exposures taken at the same 
nod position were summed after accounting for any slight spatial offsets (to within a single pixel). 

Object and telluric calibration spectra at each nod position in each echelle order were then 
extracted from a spatial profile that was > 10% of the profile peak (i.e., 3-4 pixels), using a 
background aperture off the peak (e.g., it 8-28 pixels) for residual background correction. Each 
extracted object and telluric spectrum was wavelength calibrated u sing selected telluric absorption 
lines identified from the HITRAN database (jRothman et al.lll998l ). 



Telluric features in each of the three nod positions of V1331 Cyg were removed by dividing by 
the spectrum of the telluric standard taken at the same nod position along the slit. An optimal 
telluric cancellation in V1331 Cyg was achieved using the IRAF package telluric, which allowed the 
strengths of the telluric lines to be scaled to account for a slight difference in the airmass of the 
hot standard star before normalization and division with the target spectra. 

We put our V1331 Cyg spectra on a flux scale using photometry from the literature. We 
constructed an SEP for V1331 Cyg, using the near-IR and optical photometry of 2MASS and 
Chavarrial (|l98ll ). respectively. The 2M ASS -RT-band photo metric point is in good agreement with 
the flux calibrated K-hand spectra from lNajita et al.l (|2009l ). We fit the near-IR photometry with a 
low order function to interpolate the expected flux in different L-band orders. We then rescaled our 



^See http://www2.keck.hawaii.edu/inst/nirspec/redspecl 
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NIRSPEC observations so that the flux average in each observed order matched the interpolated 
value from the SED (~ 0.35-0.41 Jy across 3.291-3.854 fim). 



3. Analysis 

The resulting L-band spectra of V1331 Cyg reveal emission features that are strong, numer- 
ous, and relatively narrow (Fig. 1). The narrow line widths are an advantage in deconstructing 
such a rich emission spectrum. The coincidence between the observed emission features and the 
wavelengths of lines of OH (red ticks) and water (blue ticks) suggest that the observed spectra are 
dominated by emission from water and OH. As we show below, the spectra are consistent with OH 
and water emission from the atmosphere of the inner disk of V1331 Cyg (Figs. 2a-h). From the 
line widths and strengths, we characterize the physical conditions in the disk where the water and 
OH emission arises. 



3.1. Disk Emission Model 

The radial velocity of the observed features has a topocentric value of -27.6 kms'^, in agree- 
ment with the radial velocity of the C^^O (2-1) emission from the envelope which provides an 



estimate of the radial velocity of the system (vlsr = — 1.3kms~ . lMcMuldrochlll993l ). The lack of 
a velocity offset with respect to the systemic velocity is consistent with an origin for the emission 
in a rotating circumstellar disk. 

We therefore modeled the observed spectra of V1331 Cyg by combining contributions from 
the stellar continuum, a disk veiling continuum, and emission lines of H2O and OH f rom a disk 



atmos phere. The disk line emission model, which follows the approach described in ICarr et al 



(|200J), assumed the geometry of a plane parallel slab of gas in a Keplerian disk that is in vertical 
hydrostatic equilibrium. The parameters describing the slab are the gas temperature (T), column 
density (A), and intrinsic line broadening, the inner and outer radii of the emission {Rin, Rout), 
and the line-of-sight disk rotational velocity at the inner radius of the emission (vsmi). 

The molecular emission spectrum can be used to constrain these parameters in the following 
schematic way. The relative strengths of the emission features constrain the temperature (T) and 
column density (A) of the emission. The large number of transitions in the spectrum and the wide 
range in their intrinsic strengths (^/-values) and energy levels are assets in this regard. The flux 
level of the emission spectrum constrains the projected emitting area for the emission, which is a 
function of Ri^ and i?outj the inclination, and the assumed distance to the system. 

The V1331 Cyg emission lines are resolved, narrow, and centrally peaked. Because the emission 
lines are resolved, we have a further constraint on the emitting radii. The maximum velocity of 
the emission profile (i.e., the velocity extent of the profile) constrains the projected velocity at the 
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inner radius of the emission (vsini). The values of Rin and f sini, along with an assumed stellar 
mass, specify the inclination of the system. 

The shape of the emission profile constrains the outer radius of the emission (i?out)) e.g., a 
small ratio of Rout/R'm produces a double-peaked profile, whereas a large ratio produces a centrally- 
peaked profile. In modeling the spectrum we find that the centrally peaked profile cannot be fit 
with Rout/Rin < 2.0; therefore the emission must arise over a range of radii and not from a narrow 
annulus. 

The temperature and column density may also vary with radius and can affect the detailed 
line profile shape. Because the line profiles are only marginally resolved in our data, our ability 
to diagnose these properties is limited. We therefore assume a constant temperature and column 
density for the emitting region in the modeling carried out here. 

Thus our strategy was to search for a combination of temperature and column density that fit 
the relative line strengths and then choose an emitting area (i.e., an outer radius) that matched 
the overall flux of the lines for the assumed distance. The combination of an observed vsini, an 
emitting area inferred from the emission flux, and a constraint on Rout/Rin from the line profile 
then allowed us to constrain the model parameters i?out! -^in! and the inclination i. We required 
that the model provide a decent fit to all four of the orders observed with the LI and L2 grating 
settings. 

The additional parameters in the modeling are the stellar mass and the distance and extinction 
to the source. F or consistency wit h our earlier study of the if-band water and CO einission 



from this source (jNajita et al.M2009l ). we adopted a distance of 550 pc (IShevchenko et al.l Il991 
Alves et al.lll998l ). a stellar mass of 1.8Mq, and no K-haiid extinction along the line of sight. 
Also consistent with our earlier study, we assumed a stellar radius and temperature of 2.2Rq 
and 7200 K, respectively, and w e estimated t h e und erlying stellar L-band continuum assuming a 



blackbody function. As noted in lNaiita et al.l (|2009l ). the stellar properties are not well known 



We constructed separate syntheses for the OH and water emission, and then combined the 
models additively. This is an acceptable approach because there is little overlap between the lines 
from the water and OH line lists to within several kms~^, despite the high density of water lines in 
the BT2 list. To determine the strength of the disk veiling continuum, we started with the stellar 
contribution to the L-band continuum and added the synthesized water emission spectrum. Disk 
continuum emission was then added to bring the model emission up to the observed level of the 
emission in each order. The required amount of disk continuum emission increased monotonically 
from ~ 0.26 Jy to ~ 0.37 Jy across L-band orders 23 to 20, consistent with the SED shape of 
V1331 Cyg (§2.2). Since most of the L-band continuum is dominated by the excess disk emission 
(i.e., veiling of approximately 10 in the L-band), no stellar spectral features would be detectable 
in our data. As the last step, a synthetic OH emission spectrum was added to produce the final 
model spectrum (Figs. 2a~h). 
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3.2. Molecular Linelists 



A successful mode l requires good line lists. Our OH line list comes from the HITRAN database 
(jRothman et al.lll998l ). an empirically determined list based on absorption in the Earth's atmo- 
sphere. Many of the OH emission lines in our L-band spectra are prominent and mostly isolated 
from one another, allowing us to identify the transitions present in our data. 



For water, we use a new theoretical line list, BT2 (jBarber et al.ll2006l ). which includes over 
500 million transitions. We trimmed the line list to a parameter range that is relevant for the 
present application, by sorting lines based on their relative strengths at 1500K, the approximate 
temperature we found in fitting the K-hand water emission of VI 331 Cyg. Here the relative strength 
is the product of the absorption cross section and the fractional population in the lower energy level 
of the transition. We excluded lines that had relative strengths < 0.001% of the strongest line in 
each order. The resulting line list included ~4600 to ~6700 lines in each of our L-band orders. 



3.3. Water and OH emission in the L-band 

We find that we can fit the water emission line profile and reproduce much of the water 
emission structure we observe with a temperature of ~ 1500 K, a line-of-sight H2O column density 
of ~ 2 X lO^^cm"^, an emitting area extending from i?in = 5.5Rq to i?out = 3.5i?in (i.e., 0.03- 
0.09AU), and vsini = 14kms~^ at Rin- The resulting emission profile is convolved with the 
3-pixel NIRSPEC slit resolution of 12.5 km s~^. We refer to the model with these parameters 
as our "reference model" for the water emission (Figs. 2a-h). As the fits demonstrate, the line 
emission is dominated by water in all observed L-band orders. The water lines detected in the 
various L-band orders span a broad range from optically thin to thick (e.g., th^o ~ 10~^ to 10^ for 
the reference model). 

These values for the water temperature and column densit y are similar to th e fit we found 



previously to the i^-band water emission observed in V1331 Cyg ([Naiita et al.ll2009l ). The emitting 
area we derive for the L-band water emission is somewhat smaller than (~ 50% of) the size of the 
K-hand emitting area found in the earlier study. This difference does not mean that the results 
are in confiict, because the K- and L-band observations were taken at different epochs and the 
line emission strength may vary. But it is also unclear whether variability is indicated because 
the present data were fiux calibrated using data from the literature rather than contemporaneous 
observations (§2). 

An inner radius of 5.5Rq is 2.5 times larger than the stellar radius of 2.2Rq (§3.1). For a 
stellar mass of 1.8Mq (§3.1), the adopted values of vsini and i?in correspond to an inclination of 
3 degrees. The inclination and the range of emission radii are consistent wit h the values we found 



in modeling the iC-band water emission from this source (iNajita et al.ll2009l ). 



The inferred values for i?in, -Rout, and inclination depend on the assumed distance to the source 
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and the stellar mass, neither of which are well known. The distance has been variously estimated 
in the range 500-800 pc (Najita et al. 2009), allowing for larger values than we assumed. A larger 
distance would require a larger disk emitting area to produce the same flux. For those larger disk 
emitting radii, a larger inclination would be needed to produce the same vsini at a fixed stellar 
mass. However, a larger distance would also affect the estimate of the stellar mass based on the 
SED, with further consequences for the Ri^, -Routi and inclination. Because of these uncertainties, 
the values of Ri^, Rout, and inclination are not well constrained by our modeling. Nevertheless, the 
modeling shows that the notional values of stellar mass and distance from the literature do allow 
for a simple disk model for the molecular emission. 

To estimate the range of column density (or temperature) that would produce a reasonable fit 
to the water emission, we held temperature (or column density) fixed at the value of the reference 
model and allowed column density (or temperature) to vary with emitting cirGSj clS cl free parameter. 
The emitting area was determined by fitting the spectrum in LI order 21, and the value was kept 
the same for all other orders. LI order 21 was chosen since it contained a mix of optically thick 
and thin lines and is a spectral region with relatively good telluric transmission. We determined 
the value of the emitting area by minimizing the sum of the ab solute values of the pixel differences 



between the observed spectrum and the model in LI order 21 (jBr anhami 1 1 982l ) . Regions with poor 
telluric transmission (i.e., < 80%) were excluded in this estimate. 



Figure 3a shows the resulting fits to LI order 21 at a fixed temperature (T=1500K) for the 
reference model (middle spectrum) and for column densities 10 times lower (upper spectrum) and 10 
times higher (lower spectrum) . The low column density model is more optically thin and has larger 
peak-to-trough variation, while the high column density model has much less contrast. The peak- 
to-trough variation in the reference model more closely resembles that of the observed spectrum. 
Figure 3b shows the model results for a neighboring order (LI order 22) obtained with the same 
emitting area used in LI order 21. The low and high column density models fit much worse than 
in LI order 21 because the emitting area was not optimized for this order. These model fits are 
clearly worse than the reference model in fitting the observed spectrum. 

Figure 4a shows fits to LI order 21 for the reference model temperature (middle spectrum) 
and for lower (upper spectrum) and higher (lower spectrum) temperatures. The lower temperature 
model has too large a peak-to-trough variation, while the higher temperature model has too little, 
compared with the reference model case (1500 K). This same trend is evident in another order (LI 
order 22, Fig. 4b). In both figures, the high and low temperature models are clearly worse than 
the reference model case. Therefore, for the water emission it appears that we can rule out column 
densities < 2 x 10^° and > 2 x 10^^ cm'^ and temperatures < 1200 K and > 2000 K. 

We adopted a similar approach in fitting the OH emission. We found that the OH emission 
lines in our L-band data can be fit with a temperature of ~ 1500 K and a line-of-sight OH column 
density of ~ 1 x lO^'^cm"^, and a larger emitting area (2.1 times larger) than the water emission, 
with radii extending from i?in = 5.5Rq to -Rout = 5i?in, and vsini = 14kms~^ at Rin- We refer to 
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the model with these parameters as our "reference model" for the OH emission (Fig. 2a-h). The 

difference between the emitting areas of water and OH in the reference model is not significant. If 
we assume the OH emission comes from the same range of radii as the water emission {Rin = 5.5Rq 
to Rout = 3.5i?in), then we find a similar fit to the OH emission at the same column density as 
above, but at a slightly higher temperature (1850 K). The OH lines in our spectra range from 
optically thin to thick, depending on the transition and model parameters (e.g., tqh ~ 0.1 to 30 
for the reference model, Table 1). The OH transitions also span a range in excitation, including 1-0 
lines with J = 11, 12, 14-18 and 2-1 lines with J = 8-18 (red ticks. Fig. 2a-h, see Table 1). Since 
the OH lines in our data span a range of optical depths and excitation levels, we can constrain the 
column density and temperature of the model fits. 

By requiring a consistent fit to all orders we can exclude column densities < 5 x 10^^ and 
> 10^^ cm~^ and temperatures < 1200 K and > 2300 K. Figure 5a compares observed spectral 
regions where OH is present (black histogram) with model fits that include H2O and OH and 
adopt higher and lower OH column densities than the reference OH model at a fixed temperature 
(T=1500K). The top panel compares the reference model (N=l x lO^^cm"^, solid blue line) with 
a low column density fit (N=5 x 10^^ cm^^: dotted red line). The low column density model 
overpredicts the flux of the lines at shorter wavelengths (3.282-3.287 ^m). The bottom panel 
compares the reference model (N=l x lO^'^cm"^, solid blue line) with a high column density fit 
(N=l X 10^^ cm~^, dotted red line). The high column density model underpredicts the flux of the 
lines at shorter wavelengths (3.282-3.287 /xm) and somewhat overpredicts the flux of the lines at 
longer wavelengths (3.776-3.780 //m). The reference model fits all the OH lines adequately well 
across a range of optical depths. 

Similarly, Figure 5b compares selected OH lines (black histogram) with fits at higher and 
lower temperatures than the reference OH model at a fixed column density (N=l x lO^'^cm"^). 
The top panel compares the reference model (T=1500 K, solid blue line) with a low temperature fit 
(T=1200K; dotted red line). The low temperature model overpredicts the flux of the low-J lines at 
shorter wavelengths (3.282-3.287 fim). The bottom panel compares the reference model (T=1500 K, 
solid blue line) with a high temperature fit (T=2300K, dotted red line). The high temperature 
model overpredicts the flux of the high-J lines at longer wavelengths (3.776-3.780 /xm). The low 
temperature model puts too much energy into the low-J lines (shorter wavelengths, top panel), and 
the high temperature model puts too much energy into the high-J lines (longer wavelengths, bottom 
panel). The reference model fits all the OH lines adequately well across a range of excitation levels. 

The good agreement between the observed spectrum and the OH-I-H2O model spectrum (Fig. 
2) shows that it is possible to account for much of the spectral structure with our reference model 
of molecular emission from a rotating disk. The temperature of the water and OH emission fit by 
the reference models is the same (~ 1500 K for both), while the OH column is ~ 10 times less than 
the water column. The range of emitting radii for both models overlaps, with the emission from 
both species arising from within 0.2 AU of the star. Some discrepancies between the observed and 
model spectra are noticeable and may be attributed to errors in the transition frequency and/or 
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oscillator strength of some lines in the BT2 water line list (e.g., features at 3.485 ^um, 3.492 ^m, 
3.499 fiui, 3.505 fim, 3.514 ^m; vertical arrows in Fig. 2c). 

The advantage of using the BT2 list, compared to the more familiar HITRAN list, when 
studying warm water emission is apparent if we compare synthetic spectra produced with the two 
line lists at a temperature of 1500 K (Fig. 6). While the fit parameters are identical in both cases, 
the BT2 list clearly matches most of the observed structure in the data, while the HITRAN list 
misses most of it. This is because of the absence of high temperature lines in HITRAN, which are 
present in the BT2 list. For water emission at lower temperatures (~ 1000 K) and in wavelength 
regions dominated by low excitation lines, the HITRAN list can do ab out as well as BT 2, as in the 
case of the 2.9 //m water emission observed in AS 205 A and DR Tau (jSalvk et al.ll2008l ). 

The line widths of individual OH and H2O lines are difficult to measure directly from our 
data owing to the crowding of numerous H2O lines throughout the L-band. In the absence of 
microturbulent broadening, the intrinsic line widths of the molecular emission we detect in our 
observations is determined by pure thermal broadening, which is ~1 kms~^ in our case for both 
H2O and OH at 1500K, and 0.8 kms"^ for CO at 1800 K. The modehng results show that the 
L-band ro-vibrational water and OH emission in this source have similar emission line profiles 
(FWHM ~ 17-20 kms~^), indicating that emission lines in both species arise from a similar range 
of disk radii. Thus it seems plausible that the two species probe similar gas temperatures. If this 
is the case, it implies that the ratio of OH and water column densities is ~ 0.1. 



3.4. Water and CO emission in the iT-band 



Water and CO emission from V1331 Cyg has been previously reported in the if -band (jNaiita et al. 



20091 ). based on an observation made in 1999, two years earlier than the observation date of the 
-L-band data reported here. To determine how the OH and water emission properties in the L-band 
compare with those of CO and water reported in the K-hand study, we re-analyzed the /C-band 
data using the same simple slab model describ ed in §3.1. Our analysis here differs in several ways 
from the previous analysis (jNajita et al.ll2009l ) in that we report a single fixed temperature and 
column density for each species rather than using radial gradients, and chemical equilibrium is not 
assumed for either species. Our analysis also uses the new BT2 line list, compared with a previous 
wate r line list that w as empirically calibrated but for a limited region near the 2-0 CO bandhead 
(i.e.,|carZeiaDl2Qo3). 



We find a good fit to the water features in the /C-band using the same parameters found for the 
L-band (i.e., ThjO = 1500 K and A^H20 = 2 x 10^^ cm^^), but with an emitting area that is 2 times 
larger (i.e., R = 5.5 — 26 Rq). It is possible that the different values we find for the emitting area 
of water in V1331 Cyg over different observational epochs result from real time variability in the 
emission, but some fraction of the difference could also be accounted for by uncertainties in the flux 
calibration in either epoch. Guided by the previous fit values for CO, we found we could fit the CO 
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emission reasonably well with a temperature of Tqo = 1800 K and a line-of-sight column density 
-^co = 6 X lO^^cm"-^, that emits over disk radii 2.3-18 Rq, with vs'm i = 23kms~^ at i?in and 
turbulent broadening of Vturb = 4km s~^ (Fig- 7). The BT2 list does a good job reproducing most 
of the emission structure observed blueward of the CO bandhead (Fig. 7), and discrepancies in the 
fit to the data may illustrate the need for empirical calibration of its line strengths and transition 
energies in the ET-band as well. The inner radius for the emission may appear uncomfortably 
close to the stellar radius {2.2Rq; §3.1). We note that the inferred disk emitting radii are not 
well constrained by our modeling and larger emitting radii are possible, in principle, with a larger 
source distance (§3.3). A more complete modeling effort that explores the uncertainty in the source 
distance is needed to explore this possibility. 

From our analysis of the iT- band data, w e find a CO-to-water column density ratio of ~ 3- 
10, consistent with the results of Najita et al. (j2009l ). The larger value would apply if the water 



emission has the same value of fturb = 4kms~^ as the CO emission. The smaller value would 
apply if the water emission has no turbulent broadening, as in the L-band modeling. The CO and 
water emission may experience different amounts of turbulent broadening if they are produced at 
different heights or radii in the disk atmosphere. We compare in Table 2 the results of the -fC-band 
modeling with the properties of the reference model for the L-band modeling. As shown in Table 
2, the properties of the water emission observed in the K- and L-bands are quite similar. 



4. Discussion 



4.1. Comparison to other sources with disk emission 

To date, gaseous wat er and OH emiss i on have been reported in the near- I R spectra of only 



a handful of young s tars (ICarr e t al. 200 4]: iNaiita et al.l l200d : iThi fc Bik 



Mandell et al.l l2008l : iNajita et al.l ,2009, : iFedele et al 
sion temperature (~ 1500 K) and column density ( 



2011 



2005 



Salyk et al. 



2008 



~ 10^^ cm-2 



As noted in §3.3, the water emis- 
that we find for V1331 Cyg are 



similar to the wate r emission properties reported for V1331 Cyg based on its /C-band spectrum 
(INaiita et al.ll2009l ). These properties are a lso similar to those found for the near-IR water emis- 
sion detected from SVS 13 (jCarr et al.ll2004l ). for which the radii of emission are comparably small 
(< 0.3 AU). The water emission detected from the disk surrounding the high mass star 08576nr292 
(6M0) was also found to have a similar temperature but a significantly lower column density 
10^^ cm~^) em itting over a larger area from an inner radius of 2 AU to an outer radius of 4 AU 



(|Thi fc Bikll200,5l ^ 



In comparison to the T Tauri stars, AS 205 and DR Tau (|Salvk et al.l 120081 ) . which also show 
water emission in the L-band, the emission spectrum of the more extreme T Tauri star V1331 Cyg, 
reported here, is much richer (i.e., has a higher density of emission lines), probably as a consequence 
of its higher water emission temperature and column density . Its spectrum is similar to that of 
SVS 13, which was previously reported by INaiita et al.l (j2007l ). 
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The Herbig Ae stars AB Aur and MWC 758 (jMandell et alj|2008) also show OH emis sion in 



the L-band, but no water emission even at a signal-to-noise of ~ 1000. iMandell et alj (|2008l ) find a 
temperature of ~ 700 K for the (optically thin) OH emission that they detect. Assuming that the 
gas is collisionally excited (and in LTE) and that the emitting region is constant in col umn density 
and extends out to ~ 15 AU, consistent with Herbig Ae disk atmosphere model of iKamp et al 



(|2005l ). the OH mass that they infer corresponds to a low OH column density of 
Lower column densities are inferred if the OH emission arises through fluorescence. 



10^^ cm" 



At Spitzer IRS wavelengths (10-40 u m) , prominent water and OH emission is commo n among T 
Tauri stars but rare among Herbig stars (jPontoppidan et al.ll2010al : ICarr &: Naiitall201ll ). Where it 
has been studied in detail, the water and OH emission detected with Spitzer is also characterized by 
larger radii and more modest temperatures and column den sities than those found for yi331 Cyg. 



T Tauri wa ter and OH emission temperatures are > 500 K (I Carr fc Najitall2008l . 



2011 



2011 



2008 . 201 ih. with column densities of ~ 10^^ cm^^ for H2O and 10^^ 



10^^ cm-2 



Carr &: Najita 



2011 



Salvk et al 



(jSalvk et al 



respective ly) for OH. The OH -to- water column density ratio of ~ 0.13 
for the typical T Tauri star AA Tau (ICarr &: Najitall2008l ) is similar to the ratio of ~ 0.1 found for 
V1331 Cyg. Thus, the temperature and, in particular, the column densities of the OH and water 
emission from V1331 Cyg are at the upper end among sources with detected water and/or OH 



emission. 



Interestingly, V1331 Cyg does not show strong water emission at Spitzer wavelegnths (ICarr &: Najita 



2011 



) indicating that the disk atmosphere at larger radii than is probed by the L-band does not 
produce much water emission. Perhaps this is because little grain settling has taken place in the 
V1331 Cyg disk atmosphere a t such radii, a situation that would limit line emission from the disk 



atmosphere (jSalvk et al 



2011 



Bethell &: BerginI ([20091) have recently pointed out that because the UV ph otoabsorption cross- 



1996 



sections of OH and water are ap proximately continuous at 10~^^ — lO^^^cm^ ( Yoshino et al 
van Dishoeck Sz Dalgarndll984l ). OH and water column densities of ~ 10^^ — 10^^ cm~^ can begin 
to shield the disk midplane from UV irradiation in the event of complete dust settling out of the 
disk atmosphere, enabling chemical synthesis to continue in the midplane despite a harsh external 
UV environment. The high column densities of OH and water that are observed in V1331 Cyg and 
other sources are more than adequate for this purpose. 



4.2. Disk Photochemistry 



^In the analysis of spectrally unresolved spectra (e.g., as observed with ^piteer IRS. ICarr &: Najita 

20081 ). the assumption is sometimes made that molecular species with similar excitation tempera- 



tures probe the same region of the disk. With that assumption, one can ratio the column density 
measured for each molecular species as a spatial average over the entire emitting region and obtain 
a rough estimate of the molecular abundances of the emitting gas. The derived abundances are 
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uncertain, in part, because it is unclear whether the emission in the different molecular species arise 
from the same region of the disk. The emission in the different species may arise at different radii 
or vertical heights in the disk atmosphere. With high resolution spectra, as employed here, we can 
probe the column density of the molecular emission as a function of disk radius for each molecular 
species, thereby providing a more robust way to probe the molecular abundances of the emitting 
gas. 

As described in section 3 , we found that th e L-band OH and water emission from V1331 Cyg 
probe a similar range of radii. ISalvk et al.l (|2008l ) found a similar result in their L-band spectroscopy 
of the T Tauri disks DR Tau and AS 205A. Ratioing the column densities found for the OH and 
water emission from these sources yields column density ratios of OH/H2O ~ 0.05-0.3. 

Thermal-chemical models of inner disk atmospheres that do not include UV irradiation typ- 
ically predict much smaller OH/H2O abundance ratios (~ 10~^. iGlassgold et al.ll2004l . |2009| ). In- 
cluding UV irradiation of the disk, which can ph otodissociate water to produce OH, can lead 
to OH/H2O ratios mor e similar to those observed (jBethell &: BerginI l2009l : IGlassgold et al.l l2009l : 
Willacy Woodsl 120091') . Perhaps UV ir radiation is also responsible for the more extreme ratio of 
OH/H2O found bv iMandell et al.l ^20od ) in their L-band spectroscopy of AB Aur and MWC 758, 
where a low column density of OH is detected in emission, but no water emission is detected. UV 
irradiation has bee n invoked to acco u nt for the water eni i ssion p roperties measured for the 6Mq 
source 08576nr292 (jThi fc Bikll2005l ). iPontoppidan et al.l (l2010al ) also discuss the possible role of 
UV photodissociation in accounting for the low equivalent width of OH and water emission (often 
consistent with non-detection) in Spitzer spectra of Herbig AeBe stars. 



In the scenario explored bv lBethell &: BerginI (|2009l ). a disk irradiated by a stronger UV field 
will enhance the column density of OH in the disk atmosphere by photodissociating water. The 
maximum OH column density is predicted to be limited to that where the OH optical depth to 
FUV photons is close to unity or an OH column density of ~ 2 x 10^'' cm~^. Since water can 
heat disk atmospheres through the absorption of FUV photons as well as provide a major coolant 



for disk atmospheres (jPontoppidar 



et al 



2010al'). it may r egulate the thermal structure of disk 



Bethell Sz Bergin (j2009l ) suggest that the water emitting layer 



atmospheres where it is abundant, 
may therefore be limited to a column density at which water is optically thick to FUV photons, 
i.e., to a water column density of < 10^*^ cm~^. These OH and water column densities are similar 



to those measured for some T Tauri disks (jCarr Sz Naiitall2008l : ISalvk et al.ll2008l ). 



The much larger column densities of both water and OH that are found for V1331 Cyg and 
SVS 13, compared to these predictions, suggest that additional processes beyond FUV heating and 
simple photochemistry play a role in the thermal-chemical properties of disk atmospheres. Since 
both V1331 Cyg and SVS 13 are high accretion rate sources, one possibility is that non-radiative 
(accretion-related) mechanical heating deepens the disk temperature inversion region and leads to 
larger column densities of molecular emission compared to more typical T Tauri stars. In models of 
disk atmospheres of typical T Tauri stars, accretion-related mechanical heating has been found to 
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enhance the column density of w arm water at the disk surface up to the level observed with Spitzer 
at AU distances from the star (jGlassgold et al.l |2009|) . The models predict much larger column 
densities of warm water < 0.25 AU than at 1 AU (Najita et al. in prep.) and may explain, in part, 
the large water column we observe. 

Another possibility is that at the higher average temperatures (and smaller disk radii, higher 
densities) from which the V1331 Cyg and SVS 13 emission originates, there are additional ionization 
or chemical processes that lead to enhanced OH and water column densities. Whether these, or 
other processes, can account for the water and OH column densities observed here is an important 
issue for the future. 
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Fig. 1. — High signal-to-noise (SNR >150), high resolution {R = A/AA = 24,000) spectra of 
V1331 Cyg spanning 4 L-band orders observed with NIRSPEC at two grating settings. Emission 
lines of water (lower ticks, blue) and OH (upper ticks, red) are evident throughout these L-band 
spectra. Lines within 10% of the strongest line in the BT2 list for each order were selected here. 
Wavelength regions where the telluric transmission was poor (i.e., < 80%) have been excised from 
the plots. 
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Fig. 2a. — A rich infrared spectrum is observed in V1331 Cyg (upper thick histogram, both panels) 
arising from emission by numerous water lines blended with OH lines (vertical red ticks, labeled). 
We model the emission across multiple L-band orders (Figs. 2a-2h), and find a good model fit 
(solid blue line) to the data by combining H2O (dashed blue line) and OH (dotted red line) model 
spectra synthesized at T = 1500 K with different column densities (2 x 10^^ and 1 x lO^'^cm"^, 
respectively) and disk emitting areas (5.5 < Rh^q/Rq < 19 5.5 < Rqh/Rq < 28, respectively). 
Telluric absorption is indicated by the spectrum of the hot star, HR 7610 (lower thick histogram, 
both panels). Wavelength regions in the data and models where the telluric transmission is below 
80% are not plotted. (A color version of the figure is available in the online journal) 
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Fig. 2b. — As in Fig. 2a but for the L2 order 21 wavelength region. 
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Fig. 2c. — As in Fig. 2a but for the L2 order 22 wavelength region. The OH line emission (red ticks, 
lower panel) is strong and optically thick, characterized hy v = 2-1 and 1-0 lines (see Table 1). 
While the combined H2O+OH emission model (solid blue) fits much of the observed structure well, 
some wavelength offsets between the reference model and the data are discernible (vertical arrows) 
along with discrepancies in the strengths of some line groups. The vertical arrows illustrate some 
examples of water lines in this order that may need empirical adjustments to their theoretically 
computed transition wavelengths and/or oscillator strengths in the BT2 line list (see §3.2). 



- 22 - 




Fig. 2d. — As in Fig. 2a but for the LI order 22 wavelength region. 
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Fig. 2e.— 



As in Fig. 2a but for the LI order 20 wavelength region. 
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Fig. 2f. — As in Fig. 2a but for the L2 order 20 wavelength region. We observe an unknown strong 
emission feature at 3.8667 /im, which could be identified as a Mg I line (Arest = 3.86748 ;um) that is 
blue shifted by 60 kms^^ relative to the disk emission we detect, perhaps as part of an outfiowing 
wind. 
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Fig. 2g. — As in Fig. 2a but for the L2 order 23 wavelength region. 
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Fig. 2h. — As in Fig. 2a but for the LI order 23 wavelength region. 
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Fig. 3a. — Model fits of water emission to the data in a portion of the LI order 21 region 
with different column density and emitting area values for a fixed temperature (T = 1500 K). 
The emitting area was chosen to optimize the fit to the model flux level in this order. The low 
column density case (A^H20=2.0 x lO^'^ cm~^, Rout = 43R0, thin red line top panel) shows a poor 
fit to the data (black histogram), having too much contrast between the strengths of the strong 
and weak lines, indicating that the lines are too optically thin. In the high column density case 
(-^H2O=2.0 X 10^^ cm~^. Rout = I2R0, thin red line lower panel), the fit to the data is also poor 
because the strong and weak lines do not show enough contrast, since they are much more optically 
thick. The center panel shows a good fit to the data (A^H20=2.0 x 10^^ cm~^. Rout = ISR©) having 
a mix of optically thin and thick lines that reproduce the strengths of both the strong and weak 
lines in this order. 
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Fig. 3b. — As in Fig. 3a but for a portion of the LI order 22 region. The high and low column 
density models fit the data more poorly here (compared to Fig. 3a) because the assumed emitting 
area is that chosen to optimize the fit to LI order 21. 
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Fig. 4a. — Model fits of water emission to the data in a portion of the LI order 21 region with 
temperature and emitting area values for a fixed column density (A^ = 2 x 10^^ cm~^). The low 
temperature case (T = 1200-R', i?out = 35-R0, thin red line top panel ) shows a poor fit to the data 
(thick black histogram), with the contrast between strengths of the strong and weak lines being 
too great. In the high temperature case (T = 2000i^, i?out = H-R©, thin red line lower panel), the 
fit to the data is also poor because the strong and weak lines do not show enough contrast. The 
center panel shows a good fit to the data (T = 1500K, Rout = 19Rq). 
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Fig. 4b.— 



As in Fig. 4a but for a portion of the L2 order 22 region. 
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Fig. 5a. — OH model fits for higher and lower column densities than the reference OH model 
at a fixed temperature (T=1500K) compared with the observed V1331 Cyg spectrum (black his- 
togram). The top panel compares the reference OH model (N=l x 10^*^ cm~^, solid blue line) with 
a low column density model (N=5 x lO^^cm"-^; dotted red line). The low column density model 
overpredicts the flux of the lines at shorter wavelengths (3.282-3.287 ^m). The bottom panel com- 
pares the reference model (N=l x lO^'^cm"^, solid blue line) with a high column density model 
(N=l X 10^^ cm~'^, dotted red line). The high column density model underpredicts the flux of the 
lines at shorter wavelengths (3.282-3.287 fim) and somewhat overpredicts the flux of the lines at 
longer wavelengths (3.776-3.780 /im). The reference model fits all the OH lines (vertical ticks, see 
Table 1) adequately well across a range of optical depths. 
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OH models with Temparature Variations 
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Fig. 5b. — OH model fits for higher and lower temperatures than the reference OH model at a 
fixed column density (N= 1 x 10^*^ cm~^) compared with the observed V1331 Cyg spectrum (black 
histogram). The top panel compares the reference OH model (T=1500K, solid blue line) with a 
low temperature model (T=1200K; dotted red line). The low temperature model overpredicts the 
flux of the low- J lines at shorter wavelengths (3.282-3.287 //m). The bottom panel compares the 
reference model (T=1500K, solid blue line) with a high temperature model (T=2300K, dotted red 
line). The high temperature model overpredicts the flux of the high- J lines at longer wavelengths 
(3.776-3.780 ^m). The low temperature model puts too much energy into the low-J lines (shorter 
wavelengths, top panel), and the high temperature model puts too much energy into the high- J 
lines (longer wavelengths, bottom panel). The reference model fits all the OH lines (vertical ticks, 
see Table 1) adequately well across a range of excitation levels. 
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Fig. 6. — A comparison of synthetic spectra at the same temperature (T = 1500 K), H2O column 
density (2 x lO^^cm"^), and emitting area (5.5 < Rh^o/Rq < 19) as was found for the reference 
model to V1331 Cyg, but generated using different water line lists. The model using the BT2 water 
list (solid blue line) fits most of the observed structure in VI 331 Cyg (black histogram), while the 
model from the HITRAN water line list (dotted red line) does a poor job fitting the data (thick 
black histogram) , with very few emission features present arising from a depressed continuum level. 
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Fig. 7. — A fit to the i^T-band spectrum of V1331 Cyg observed in an earlier epoch (thick histogram, 
Najita et al. 2009), showing strong water emission blueward of the overtone CO emission bandhead. 
The model fit, a combination of water and CO emission (thin blue line), does a good job fitting the 
observed structure, highlighting the utility of the BT2 list in fitting warm K-hand water emission. 
The fit parameters for the water emission in the X-band are the same as for the water in the L-band 
(i.e.. Fig. 2), except for having an emitting area that is 1.6 times larger, which could be ascribed to 
time variable emission and/or uncertainties in the flux calibration between the two observational 
epochs. The CO model has a temperature of 1800 K, a column density of 6 x 10^^ cm~^, and an 
emitting area extending from Ri^ = 2.3Rq to -Rout = 8i?in (see Table 2). 
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Table 1. Physical properties of OH lines detected in our observations 



L-band Order 


Transition A (A) 


Line Transition Details 


Model Line Opacity'^ 


LI order 23 


32823.671 


i/=l-0; J=11.5e 


22.9 




32846.868 


i/=l-0; J=12.5e 


26.0 




32859.990 


j/=l-0; J=12.5f 


25.8 




32872.299 


i/=2-l; J=8.5f 


2.5 




32878.437 


^=2-1; J=8.5e 


2.5 




32926.693 


i/=2-l; J=9.5e 


2.9 


L2 order 23 


33367.886 


u=l-0; J=12.5f 


17.2 




33378.685 


i/=l-0; J=12.5e 


17.0 




33441.684 


u=2-l; J=10.5f 


2.5 




33651.661 


iy=3-2; J=7.5e 


0.16 


LI order22 


34484.234 


u=2-l; J=11.5f 


1.4 




34494.249 


^=2-1; J=11.5e 


1.4 




34573.554 


iy=l-0; J=14.5e 


8.4 




34584.036 


i>=l-0; J=15.5e 


9.3 




34602.034 


i/=l-0; J=15.5f 


9.2 


L2 order 22 


35077.159 


j/=2-l; J=12.5f 


1.1 




35088.615 


^=2-1; J=12.5e 


1.1 




35110.880 


i/=2-l; J=13.5e 


1.2 




35126.505 


u=2-l; J=13.5f 


1.2 




35132.676 


iy=3-2; J= 9.5f 


0.11 




35201.132 


i/=l-0; J=15.5f 


5.6 




35216.445 


i^=l-0; J=15.5e 


5.5 




35223.349 


u=l-0; J=16.5e 


6.1 




35243.170 


u=l-0; J=16.5f 


6.0 


LI order 21 


35874.842 


u=l-0; J=16.5f 


3.6 




35891.844 


iy=l-0; J=16.5e 


3.5 




35895.358 


i/=l-0; J=17.5e 


3.9 




35917.113 


i/=l-0; J=17.5f 


3.8 


L2 order 21 


36582.769 


u=l-0; J=17.5f 


2.2 




36601.565 


i/=l-0; J=17.5e 


2.2 




36601.815 


iy=l-0; J=18.5e 


2.4 




36625.624 


i/=l-0; J=18.5f 


2.4 


LI order 20 


37788.009 


u=2-l; J=16.5e 


0.24 




37793.784 


u=2-l; J=17.5e 


0.27 




38109.260 


u=l-0; J=19.5f 


0.75 


L2 order 20 


38532.612 


u=2-l; J=17.5f 


0.16 




38552.724 


iy=2-l; J=17.5e 


0.15 




38554.918 


^=2-1; J=18.5e 


0.17 




38580.509 


u=2-l; J=18.5f 


0.17 



a OH reference model: T = 1500 K, AT = 1.0 x lO^" cm-2 
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Table 2. Comparison of molecular gas parameters derived for V1331 Cyg 



Obs. Date 


A 


Molecule 


T 


N 


Emitting Radii 


I'turb 


A''moi/-^'^H2 


(UT) 






(K) 


(xl02lcm-2) 


(-R©) 


(kms-1) 




1999 July 3 


K 


CO 


1800 


6 


2.3 - 18 


4 


3-10 


1999 July 3 


K 


H2O 


1500 


0.6 


5.5 - 26 


4 


1 


1999 July 3 


K 


H2O 


1500 


2 


5.5 - 26 





1 


2001 July 11 


L 


H2O 


1500 


2 


5.5 - 19 





1 


2001 July 11 


L 


OH 


1500 


0.1 


5.5 - 28 





0.05 



